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Abstract 

We consider finite temperature SU (2) gauge theory in the continuum formulation, which necessi¬ 
tates the choice of a gauge fixing. Choosing the Landau gauge, the existing gauge copies are taken into 
account by means of the Gribov-Zwanziger (GZ) quantization scheme, which entails the introduc¬ 
tion of a dynamical mass scale (Gribov mass) directly influencing the Green functions of the theory. 
Here, we determine simultaneously the Polyakov loop (vacuum expectation value) and Gribov mass 
in terms of temperature, by minimizing the vacuum energy w.r.t. the Polyakov loop parameter and 
solving the Gribov gap equation. Inspired by the Casimir energy-style of computation, we illustrate 
the usage of Zeta function regularization in finite temperature calculations. Our main result is that 
the Gribov mass directly feels the deconfinement transition, visible from a cusp occurring at the same 
temperature where the Polyakov loop becomes nonzero. In this exploratory work we mainly restrict 
ourselves to the original Gribov-Zwanziger quantization procedure in order to illustrate the approach 
and the potential direct link between the vacuum structure of the theory (dynamical mass scales) and 
(de)confinement. We also present a first look at the critical temperature obtained from the Refined 
Gribov-Zwanziger approach. Finally, a particular problem for the pressure at low temperatures is 
reported. 
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1 Introduction 


Within SU{N) Yang-Mills gauge theories, it is well accepted that the asymptotic particle spectrum does 
not contain the elementary excitations of quarks and gluons. These color charged objects are confined 
into color neutral bound states: this is the so-called color confinement phenomenon. It is widely believed 
that confinement arises due to non-perturbative infrared effects. Many criteria for confinement have been 
proposed (see the nice pedagogical introduction m)- A very natural observation is that gluons (due to 
the fact that they are not observed experimentally) should not belong to the physical spectrum in a 
confining theory. On the other hand, the perturbative gluon propagator satisfies the criterion to belong 
to the physical spectrum (namely, it has a Kallen-Lehmann spectral representation with positive spectral 
density). Hence, non-perturbative effects must dress the perturbative propagator in such a way that 
the positivity conditions are violated, such that it does not belong to the physical spectrum anymore. 
A well-known criterion is related to the fact that the Polyakov loop [5] is an order parameter for the 
confinement/deconfinement phase transition via its connection to the free energy of a (very heavy) quark. 
The importance to clarify the interplay between these two different points of view (nonperturbative Green 
function’s behaviour vs. Polyakov loop) can be understood by observing that while there are, in principle, 
infinitely many different ways to write down a gluon propagator which violates the positivity conditions, 
it is very likely that only few of these ways turns out to be compatible with the Polyakov criterion. 

One of the most fascinating non-perturbative infrared effects is related to the appearance of Gribov 
copies [3] which represent an intrinsic overcounting of the gauge-field configurations which the perturba¬ 
tive gauge-fixing procedure is unable to take care of. Soon after Gribov’s seminal paper. Singer showed 
that any true gauge condition, as the Landau gaug^ presents this obstruction |3] (see also 0). The 
presence of Gribov copies close to the identity induces the existence of non-trivial zero modes of the 
Faddeev-Popov operator, which make the path integral ill defined. Even when perturbation theory 
around the vacuum is not affected by Gribov copies close to the identity (in particular, when YM-theory 
is defined over a flat space-tim^ with trivial topology m), Gribov copies have to be taken into ac¬ 
count when considering more general cases (such as with toroidal boundary conditions on flat space-time 
[mill]). Thus, in the following only the standard boundary conditions will be considered. 

The most effective method to eliminate Gribov copies, at leading order proposed by Gribov himself, 
and refined later on by Zwanziger m usi na m] corresponds to restricting the path integral to the 
so-called Gribov region, which is the region in the functional space of gauge potentials over which the 
Faddeev-Popov operator is positive definite. The Faddeev-Popov operator is Hermitian in the Landau 
gauge, so it makes sense to discuss its sign. In [isiin] Dell’Antonio and Zwanziger showed that all the 
orbits of the theory intersect the Gribov region, indicating that no physical information is lost when 
implementing this restriction. Even though this region still contains copies which are not close to the 
identity [18], this restriction has remarkable effects. In fact, due to the presence of a dynamical (Gribov) 
mass scale, the gluon propagator is suppressed while the ghost propagator is enhanced in the infrared. 
More general, an approach in which the gluon propagator is “dressed” by non-perturbative corrections 
which push the gluon out of the physical spectrum leads to propagators and glueball masses in agreement 
with the lattice data na HO). With the same approach, one can also solve the old problem of the 
Gasimir energy in the MIT-bag model [21]. Moreover, the extension of the Gribov gap equation at 
finite temperature provides one with a good qualitative understanding, already within the semiclassical 
approximation, of the deconfinement temperature as well as of a possible intermediate phase in which 
features of the confining phase coexist with features of the fully deconfined phase in agreement with 
different approaches (see |22] and references therein). Furthermore, within this framework the presence 
of the Higgs field |23l [24] as well of a Ghern-Simons term in 2-1-1 dimensions [25] can be accounted for as 
well. 

For all these reasons, it makes sense to compute the vacuum expectation value of the Polyakov loop 
when we eliminate the Gribov copies using the Gribov-Zwanziger (GZ) approach. Related computations 
are available using different techniques to cope with nonperturbative propagators at finite temperature, 
see e.g. [MinziiMi Hi isnjisijiniiiMiisiiss ES] • In the present paper, we will perform for the first time (to 

^We shall work exclusively with the Landau gauge here. 

^In the curved case, the pattern of appearance of Gribov copies can be considerably more complicated: see in particular 
[ilTlElIl]. Therefore, only the flat case will be considered here. 
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the best of the authors’ knowledge) this computation, using two different techniques, to the leading one- 
loop approximation. In [sziiMiisn], it was already pointed out that the Gribov-Zwanziger quantization 
offers an interesting way to illuminate some of the typical infrared problems for finite temperature gauge 
theories. 

In Section 2, we provide a brief technical overview of the Gribov-Zwanziger quantization process and 
eventual effective action. In the following Section 3, the Polyakov loop is introduced into the GZ theory 
via the background field method, building on work of other people [m nil Eg. Next, Section 4 handles 
the technical computation of the leading order finite temperature effective action, while in Section 5 we 
discuss the gap equations, leading to our estimates for both Polyakov loop and Gribov mass. The key 
finding is a deconfinement phase transition at the same temperature at which the Gribov mass develops a 
cuspy behaviour. We subsequently also discuss the pressure and energy anomaly. Due to a problem with 
the pressure in the GZ formalism (regions of negativity), we take a preliminary look at the situation upon 
invoking the more recently developed Refined Gribov-Zwanziger approach. We summarize in Section 7. 


2 A brief summary of the Gribov-Zwanziger action in Yang- 
Mills theories 


Let us start by giving a short overview of the Gribov-Zwanziger framework [3i [m EH [ig. As already 
mentioned in the Introduction, the Gribov-Zwanziger action arises from the restriction of the domain of 
integration in the Euclidean functional integral to the Gribov region fl, which is defined as the set of 
all gauge field configurations fulfilling the Landau gauge, = 0, and for which the Faddeev-Popov 

operator is strictly positive, namely 

L! = {A“ ; = 0 ; - gr^'^A^) > 0 } . 

The boundary dVi of the region D is the (first) Gribov horizon. 

One starts with the Faddeev-Popov action in the Landau gauge 

5'fp = S'ym + ‘S'gf , (1) 

where S'ym and Sgf denote, respectively, the Yang-Mills and the gauge-fixing terms, namely 

Sym =IJ d^x , ( 2 ) 

and 

Sgf = J , (3) 

where (c'*, c“) stand for the Faddeev-Popov ghosts, 6“ is the Lagrange multiplier implementing the Landau 
gauge, - gP^^A^) is the covariant derivative in the adjoint representation of SU{N), and 

denotes the field strength: 

f;, = d,Al - d,Al + gr^^A^Al . (4) 

Following El EH El Eg, the restriction of the domain of integration in the path integral is achieved by 
adding to the Faddeev-Popov action Spp an additional term H{A), called the horizon term, given by the 
following non-local expression 

H{A,^)=g^J d^xd^yr^^Alix) [M-H-f)Y\x,y)f^^Al{y) , (5) 

where stands for the inverse of the Faddeev-Popov operator. The partition function can then be 

written as HElElEg: 

Zgz = [ VAVcVcVbe~^^^ = / X>c X>c , (6) 

Jn J 


3 


where V is the Euclidean space-time volume. The parameter 7 has the dimension of a mass and is known 
as the Gribov parameter. It is not a free parameter of the theory. It is a dynamical quantity, being 
determined in a self-consistent way through a gap equation called the horizon condition [3 m m ng, 
given by 

=4V(N^-1) , (7) 

where the notation {H{A, 7 ))q 2 means that the vacuum expectation value of the horizon function H{A, 7 ) 
has to be evaluated with the measure defined in Eq.(|^. An equivalent all-order proof of eq.(|^ can be 
given within the original Gribov no-pole condition framework [3], by looking at the exact ghost propagator 
in an external gauge field m- 

Although the horizon term H{A,^), eq.(|^, is non-local, it can be cast in local form by means of the 
introduction of a set of auxiliary fields ( 0 )“^, where (plf) are a pair of Bosonic fields, 

while are anti-commuting. It is not difficult to show that the partition function Zqz in eq.(|^ 

can be rewritten as US El [15] 

^GZ = , ( 8 ) 

where $ accounts for the quantizing fields. A, c, c, b, cu, oj, (p, and (p, while S’czl'h] is the Yang-Mills 
action plus gauge fixing and Gribov-Zwanziger terms, in its localized version. 


with 


and 


Sgz = S'ym + -S'gf + So + Sj , 

So = J d^x X + 


S^= 7^1 d^x + vll) - 47"E(iV2 - 1) ■ 

It can be seen from ([^ that the horizon condition ([^ takes the simpler form 

dS„ 


97 " 


= 0 


which is called the gap equation. The quantity £^( 7 ) is the vacuum energy defined by 

= Zgz . 


(9) 

( 10 ) 

( 11 ) 


( 12 ) 


(13) 


The local action ^qz in eq. is known as the Gribov-Zwanziger action. Remarkably, it has been 
shown to be renormalizable to all orders [T3E1IIS1I1I1I111I13I11I1SI- This important property of 
the Gribov-Zwanziger action is a consequence of an extenstive set of Ward identities constraining the 
quantum corrections in general and possible divergences in particular. In fact, introducing the nilpotent 
BRST transformations 



= -Dfc^ 



SC^ 


sc^ 

II 

- ab 

—ab 


= T’m > 

ab 

ab 

SV’n 

= w , 


c c 


= 0 , 
= 0 , 
= 0 , 


= 0 


50; 


(14) 


it can immediately be checked that the Gribov-Zwanziger action exhibits a soft breaking of the BRST 
symmetry, as summarized by the equation 


s-S'gz = 7^^ , 


where 


A = 


/ 74_ / rahc f TAQ-'fyi 7n\/ be . —bc\ . rabc ai 

d X [-gf [D^ c )(vJ^ +^^,)+gf A 


abc A a, ,bc\ 


(15) 

(16) 
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Notice that the breaking term A is of dimension two in the fields. As such, it is a soft breaking and the 
ultraviolet divergences can be controlled at the quantum level. The properties of the soft breaking of the 
BRST symmetry of the Gribov-Zwanziger theory and its relation with confinement have been object of 
intensive investigation in recent years, see [i51IT7lli51li^l5Ull5n[5^[551[Ml[55| . Here, it suffices to mention 
that the broken identity (151 is connected with the restriction to the Gribov region H. However, a set 


of BRST invariant composite operators whose correlation functions exhibit the Kallen-Lehmann spectral 
representation with positive spectral densities can be consistently introduced [55]. These correlation 
functions can be employed to obtain mass estimates on the spectrum of the glueballs muoi. 

Let us conclude this brief review of the Gribov-Zwanziger action by noticing that the terms Sg{ and 
So in expression © can be rewritten in the form of a pure BRST variation, i.e. 


S', 


Sn = 


* J d^x , 


(17) 


so that 


Sgz = Sym + s 








(18) 


from which eq.(15l becomes apparent. 


3 The Polyakov loop and the background field formalism 


In this section we shall investigate the confinement/deconfinement phase transition of the SU{2) gauge 
field theory in the presence of two static sources of (heavy) quarks. The standard way to achieve this 
goal is by probing the Polyakov loop order parameter, 

V = tr dt Ao(t,x)^ ^ (^ 9 ) 

with P denoting path ordering, needed in the non-Abelian case to ensure the gauge invariance of P. 
This path ordering is not relevant at one-loop order, which will considerably simplify the computations 
of the current work. In analytical studies of the phase transition involving the Polyakov loop, one usually 
imposes the so-called “Polyakov gauge” on the gauge field, in which case the time-component Aq becomes 
diagonal and independent of (imaginary) time. This means that the gauge field belongs to the Gartan 
subalgebra. More details on Polyakov gauge can be found in |2H1 EZl EH] • Besides the trivial simplification 
of the Polyakov loop, when imposing the Polyakov gauge it turns out that the quantity (Ag) becomes a 
good alternative choice for the order parameter instead of P. This extra benefit can be proven by means 
of Jensen’s inequality for convex functions and is carefully explained in |28j . see also [271 Hi 1301 ED El]- 
For example, for the SU{2) case we have the following: if ^gP {Aq) = ^ then we are in the “unbroken 
symmetry phase” (confined or disordered phase), equivalent to {P) = 0; otherwise, if \gP{Ao) < 
we are in the “broken symmetry phase” (deconfined or ordered phase), equivalent to {P) 7 ^ 0. Since 
P oc e~^'^ with T the temperature and F the free energy of a heavy quark, it is clear that in the 
confinement phase, an infinite amount of energy would be required to actually get a free quark. The 
broken/restored symmetry referred to is the Zjv center symmetry of a pure gauge theory (no dynamical 
matter in the fundamental representation). 

A slightly alternative approach to access the Polyakov loop was worked out in [31] . In order to probe 
the phase transition in a quantized non-Abelian gauge field theory, we use, following |32j . the Background 
Field Gauge (BFG) formalism, detailed in general in e.g. [5S]. Within this framework, the effective gauge 
field will be defined as the sum of a classical field and a quantum field A^: a^(a;) = a“ (a;)t“ = A^-|-A^, 
with F the infinitesimal generators of the SU{N) symmetry group. The BFG method is a convenient 
approach, since the tracking of breaking/restoration of the Zjv symmetry becomes easier by choosing the 
Polyakov gauge for the background field. 

Within this framework, it is convenient to define the gauge condition for the quantum field, 

D^Af, = 0 , ( 20 ) 
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known as the Landau-DeWitt (LDW) gauge fixing condition, where is the 

background covariant derivative. After integrating out the (gauge fixing) auxiliary field 6“, we end up 
with the following Yang-Mills action, 


5'bfg 



_ Tpo, pa 
4 


{DAf 

2 ^ 


+ eDfD^ 



( 21 ) 


Notice that, concerning the quantum field A^, the condition (201 is equivalent to the Landau gauge, 
yet the action still has background center symmetry. The LDW gauge is actually recovered in the limit 
^ > 0, taken at the very end of each computation. 


It is perhaps important here to stress that we are restricting our analysis to the (background) Landau 
gauge, for which a derivation argument in favor of the action (211 can be provided. For a vanishing 
background, this is precisely the original Gribov-Zwanziger construction mu [13, also applicable to 
the Coulomb gauge. More recently, it was also generalized to the SU(2) maximal Abelian gauge in [55] - 
Intuitively, it might be clear that the precise influence on the quantum dynamics by Gribov copies can 
strongly depend on the chosen background, given that Gribov copies are defined via the zero modes of 
the Faddeev-Popov operator of the chosen gauge condition, which itself explicitly depends on the chosen 
background. This is open to further research, as it has not been pursued in the literature yet. Though, 
for a constant background as relevant for the current purposes, it will be discussed elsewhere that the 
action is indeed obtainable via a suitable extension of the arguments of mu El]. 


In the absence of a background, a proposal for a generalization to the linear covariant gauges was 
put forward in misi], albeit leading to a very complicated nonlocal Lagrangian structure, containing 
e.g. reciprocals and exponentials of fields. To our knowledge, no practical computations were done so 
far with this formalism. Nonetheless, potential problems with gauge parameter dependence of physical 
quantities were discussed in misi], not surprisingly linked to the softly broken BRST symmetry, see 
also our Section 2 for more on this and relevant references. 


A very recent alternative for the linear covariant gauges was worked out in |62] . partially building on 
earlier work of [53]. With this proposal, it was explicitly checked at one loop that the Gribov parameter 
7 ^ and vacuum energy are gauge parameter independent. This at least suggests that in this class of 
covariant gauges, an approach to Gribov copies can be worked out that is compatible with gauge parameter 
independence |Mj. 

As explained for the simple Landau gauge in the previous section, the Landau background gauge 
condition is also plagued by Gribov ambiguities, and the Gribov-Zwanziger procedure is applicable also 
in this instance. The starting point of our analysis is, therefore, the GZ action modified for the BFG 
framework (see |M]): 


“^GZ+PLoop — I d X 




D“(aX - " 1)} • (22) 


As mentioned before, with the Polyakov gauge imposed to the background field A^, the time-component 
becomes diagonal and time-independent. In other words, we have A^(x) = Aq^^o, with Aq belonging to 
the Gartan subalgebra of the gauge group. For instance, in the Gartan subalgebra of SU{2) only the 
generator is present, so that Ag = (5 “^Aq = ^“^Ao. As explained in [32], at leading order we then simply 
find, using the properties of the Pauli matrices, 

V = cos^, (23) 

where we defined 

r = g/SAo , (24) 

with /3 the inverse temperature. Just like before, r = tt corresponds to the confinement phase, while 
0 < r < TT corresponds to deconfinement. With a slight abuse of language, we will refer to the quantity 
r as the Polyakov loop hereafter. 
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Since the scope of this work is limited to one-loop order, only terms quadratic in the quantum fields in 
the action (221 shall be considered. One then immediately gets an action that can be split in term coming 
from the two color sectors: the 3rd color direction, called Cartan direction, which does not depend on 
the parameter r; and one coming from the 2x2 block given by the 1st and 2nd color directions. This 
second 2x2 color sector is orthogonal to the Cartan direction and does depend on r. The scenario can 
then be seen as a system where the vector field has an imaginary chemical potential irT and has isospins 
-fl and —1 related to the 2x2 color sector and one isospin 0 related to the 1x1 color sector. 


4 The finite temperature effective action at leading order 


Considering only the quadratic terms of (22), the integration of the partition function gives us the 


following vacuum energy at one-loop order, defined according to (13), 


= - 


d{N^ - 1 ) 

2Ng^ 


A"^ -f ;^((i — 1) trln 


D^ + X* 

-D2 


- 2 , 


(25) 


where V is now just the spacial volume. Here, V is the covariant derivative in the adjoint representation 
in the presence of the background Hq field and = 2Ng^j'^. Throughout this work, it is always tacitly 
assumed we are working with N = 2 colors, although we will frequently continue to explicitly write N 
dependence for generality. Using the usual Matsubara formalism, we have that 'D'^ = {2'KnT -frsT)^ -f^, 
where n is the Matsubara mode, g is the spacelike momentum component, and s is the isospin, given by 
— 1, 0, or -fl for the SU{2) cas^ 

The general trace is of the form 


1 

w 


+ O0 

trln(—-f 

s n——oo 


^"2 ((27rnT -f rsTf + f + m^) , 


(26) 


which will be computed immediately below. 


4.1 The sum-integral: 2 different computations 


We want to compute the following expression: 

+00 , ^ 3 - 


^ = ^ E J j^Jr,{(2nnT + rTr+g^+m^) . 


(27) 


One way to proceed is to start by deriving the previous expression with respect to m^. Then, one can 
use the well-known formula from complex analysis 


-|-oo 


E E cot(7rz)/(2 


(28) 


where the sum is over the poles Zq of the function f{z). Subsequently we integrate with respect to m? 
(and determine the integration constant by matching the result with the known T = 0 case). Finally one 
can split off the analogous T = 0 trace (which does not depend on the background field) to find 


1 = 




( 27 r) 


^ ln(g^ -f m^) + T 


d?q 

(^ 


In 1 -f e 


-2 




- 2e 




cosr 


(29) 


where the limit e —> 0 was taken in the (convergent) second integral. The first term in the r.h.s. is the 
(divergent) zero temperature contribution. 


^The SC/(3) case was handled in m as well (see also EZl). 
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Another way to compute the above integral is by making use of Zeta function regularization techniques, 
which are particularly useful in the computation of the Casimir energy in various configurations see 
[551IB5] . The advantage of this second technique is that, although it is less direct, it provides one with an 
easy way to analyze the high and low temperature limits as well as the small mass limit, as we will now 
show. Moreover, within this framework, the regularization procedures are often quite transparent. One 
starts by writing the logarithm as In a; = — lims_>o dsX~^, after which the integral over the momenta can 
be performed: 


I = —T lim ds 

s—>-0 



r(s-3/2) 

87r5r(s) 


[( 27 rnT + rTf + ' 


(30) 


where the renormalization scale ^ has been introduced to get dimensional agreement for s ^ 0, and where 
we already put e = 0, as s will function as a regulator — i.e. we assume s > 3/2 and analytically continuate 
to bring s —>■ 0. Using the integral representation of the Gamma function, the previous expression can 
be recast to 


X = —T lim ds I —3- / 


^s-5/2g-t((27rnT+rT)^+m^)^^ 


= — lim ds n 


2s 


n——oo 
rn4—2s 


s—>-0 ^ I 4s^2s 


dyy' 




~e e 

n——oo 


-V{n+^) 


(31) 


where the variable of integration was transformed as y = > 0 in the second line. Using the 

Poisson rule (valid for positive ui): 


+ O0 

E 


^-{n+xfui _ 


n— — oo 


i + ^E' 


cos {2mrx) 


(32) 


n—1 


we obtain that 


X = — lim dsU^ 

s-s-O 


’r(s - 2)T4-2® 

/ \ 

4®7r2®-2r(s) 

iv47r2T2 ) 


2-s 


+ 


P'4 


/ 9 \ 1 — 3/2 OO 


4® '^7r®r(s) \4Tr‘^T‘^ J 


n—1 


\ T J 


where K^{z) is the modified Bessel function of the second kind. Simplifying this, we find 


X = 


2(47r)2 


f Tn^\ 

3' 

In ( —rr 


1. vmv 

~ 2 


-E 


n—1 


m?T^ cos {nr) /nm\ 

)/2/)2 ^2 [-7^) , 


, (33) 


(34) 


where the first term is the T = 0 contribution, and the sum is the finite-temperature correction. Using 
numerical integration and series summation, it can be checked that both results (291 and (341 are indeed 
identical. Throughout this paper, we will mostly base ourselves on the expression (29). Nonetheless the 
Bessel series is quite useful in obtaining the limit cases m = 0, T —)■ oo, and T —)■ 0 by means of the 
corresponding behaviour of K 2 {z). Observing that 


lim — 

m—^0 \ 


m^T‘^K 2 (^) cos(nrs) \ 2T^ cos(nrs) 


(35) 


we obtain 


p4 


I™=0 = - ^ [Li4 + LU , 

where Li 3 (z) = is the polylogarithm or Jonquiere’s function. 

Analogously, 


(36) 


lim K 2 ( 


fmn\ 


2T^ 


1 


T^oo ^ \ T J 2 ’ 


























so that 


^T—^oo — 


m 


( w?\ 

3' 

In ^ 


1. VMV 

~ 2 


+ ^ {m^ [Li2 + Li2 (e*"^)] - [m (e-*"^) + m (e*"*)] } . 


2(47r)2 

Finally for T —>■ 0 we can use the asymptotic expansion of the Bessel function ITS]: 

(f, 


ak{v)\ , . , M 3 

, |Arg(^)|<-^, 


\k^Q 


where afc(z^) are finite factors. So, at first order (A: = 0), 


^T-)-0 — 


m 






2(47r)2 [ V / 2j 2^2713/2 L = 

4.2 The result for further usage 


Lis (e-T-i’-«) +Li 5 


Making use of the result (291 we may define 


I{m'^,r,s,T) =T / In f 1 + e 2^'^"^’" _ 2e cosrs 


(27r)s 


so that the vacuum energy (251 can be rewritten as 

£v = - 1)(^^ - l)trr=o In ^^g2 ~ ~ 1) ti'r=o ln(-9^) 


2Ng^ 


+ X! +/(-*A^,r,s,r) - 7(0, r,s,T)) - i/(0,r,s,r)^ 


(37) 

(38) 

(39) 


(40) 


(41) 


where tr^^o denotes the trace taken at zero temperature. 


5 Minimization of the effective action, the Polyakov loop and 
the Gribov mass 

5.1 Warming-up exercise: assuming a T-independent Gribov mass A 


As a first simpler case, let us simplify matters slightly by assuming that the temperature does not 
influence the Gribov parameter A. This means that A will be supposed to assume its zero-temperature 
value, which we will call Aq, given by the solution of the gap equation 0 at zero temperature. In this 
case, only the terms with the function I matter in (411, since the other terms do not explicitely depend 
on the Polyakov line r. Plotting this part of the potential (see Figure]^, one finds by visual inspection 
that a second-order phase transition occurs from the minimum with r = tt to a minimum with r 7 ^ tt. 
The transition can be identified by the condition 


Using the fact that 


2 


dr2 


,T) = -2T 


= 0 


d^q 




(27r)3 


1 -I- e’ 


. v5 


(42) 


(43) 


when s = ±1 and zero when s = 0, the equation (421 can be straightforwardly solved numerically for the 
critical temperature. We find 

rcrit=0.45Ao. (44) 
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[£„(T)-£„(T=0)]/aJ 



Figure 1: The effective potential (41) at the temperatures (from below upwards at r = tt) 0.42, 0.44, 0.46, 
and 0.48 times A as a function of r, with the simplifying assumption that A maintains its zero-temperature 
value Ao throughout. It can be seen that the minimum of the potential moves away from r = tt in between 
T = 0.44 A and 0.46 A. 


5.2 The T-dependence of the Gribov mass A 

Let us now investigate what happens to the Gribov parameter A when the temperature is nonzero. Taking 


the derivative of the effective potential (411 with respect to A^ and dividing by d{N^ — 1)\^/Ng^ (as we 


are not interested in the solution A^ = 0) yields the gap equation for general number of colors N-. 

1 Id-I Ng^ 


Id- 1 2 

~2d ^ ^’’a4-hA4 ' 2 d iV2-lA2 




( ^2 


,T) 


(45) 


where the notation dl / dm^ denotes the derivative of / with respect to its first argument (written in 
(401). If we now define Ag to be the solution to the gap equation at T = 0: 


, Id - 1 2 1 


(46) 


then we can subtract this equation from the general gap equation (451. After dividing through (d — 
VfNg"^/2d and setting d = 4 and N = 2, the result is 

E T)]=0, (47) 


1 


1 


d'^q 

+ ~ q^ + Xt 


dm? 


where now all integrations are convergent. This equation can be easily solved numerically to yield A as 
a function of temperature T and background r, in units Ag. This is shown in Figure 


5.3 Absolute minimum of the effective action 

As A does not change much when including its dependence on temperature and background, the transition 


is still second order and its temperature is, therefore, still given by the condition (42l. Now, however. 


the potential depends explicitely on r, but also implicitely due to the presence of the r-dependent A. We 
therefore have 


dr2 " 


d'^Sv ^ ^dX d'^Ey 
F 2 — — —— “1“ 


9r2 


dr drdX 


d'^X d£y 

fdX\ 

' d^Sy 

dr'^ dX 

\drj 

9A2 


(48) 


X=\{r),r=Tr 


Now, dA/dr|r=7r = 0 due to the symmetry at that point. Furthermore, as we are considering A ^ 0, 
dSy/dX = 0 is the gap equation and is solved by A(r). Therefore, we find for the condition of the 
transition: 

d'^S 

^(r,A,r) =0, (49) 
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T/A,, 


Figure 2: The Gribov parameter A as a function of the temperature T at r equals to zero (upper line) 
and TT (lower line), in units of the zero-temperature Gribov parameter Ag. 


where the derivative is taken with respect to the explicit r only. 


We already solved equation (491 in section |5.1[ giving (|44[) : 


T = 0.45A(r,T) . 


(50) 


As we computed A as a function of r and T in section |5.2| already, it is again straightforward to solve this 
equation to give the eventual critical temperature: 


Fcrit — 0.40 Aq , 


as expected only slightly different from the simplified estimate (441 found before. 


(51) 


5.4 The T-dependence of the Polyakov loop r and the equation of state 

5.4.1 Deconfinement transition and its imprint on the Gribov mass 


Let us now investigate the temperature dependence of r. The physical value of the background field r is 
found by minimizing the vacuum energy: 


-—£y = 0 . 

ar 


From the vacuum energy (411 we have 
d£, 


dr 


= {d-l) 


■ 2 , dl 2 d dl 


(0,r,r) 


= 0 . 


(52) 


(53) 


The expression (531 was obtained after summation over the possible values of s. Furthermore, we used 
the fact that I{rn?^r, 4-1, T) = I{m?, r, —1, T) and that s = 0 accounts for terms independent of r, which 
are cancelled by the derivation w.r.t. r. From (401 one can get, whenever s = ±1: 


dl{rn^, r, T) 
dr 


= T 


d^q 


2e 




sinr 


(27r)3 


1 + e 


-2 




- 2e 




cosr 


(54) 


Since ( |53[ ) is finite, we can numerically obtain r as a function of temperature. From the dotted curve in 
Figure l^one can easily see that, for T > Tcrit ~ 0.40Ao, we have r ^ n, pointing to a deconfined phase, 
confirming the computations of the previous section. In the same figure, A(T) is plotted in a continuous 
line. We observe very clearly that the Gribov mass A(T) develops a cusp-like behaviour exactly at the 
critical temperature T = Tcrit- 
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Figure 3: The dotted line curve represents r(T), while the continuous line is A(T). At T « 0.40Ao, both 
curves clearly have a discontinuous derivative. 

5.4.2 Equation of state 


Following m, we can also extract an estimate for the (density) pressure p and the interaction measure 
I/T*, shown in Figure]^ (left and right respectively). As usual the (density) pressure is defined as 


P = 


1 

w 


In^Gz , 


(55) 


which is related to the free energy by p = —Here the plot of the pressure is given relative to the Stefan- 
Boltzmann limit pressure: psB = where k = {N'^ — l)7rT^/45 is the Stefan-Boltzmann constant 
accounting for all degrees of freedom of the system at high temperature. We subtract the zero-temperature 
value, such that the pressure becomes zero at zero temperature: p{T) = — [f„(T) — £„(T = 0)]. Namely, 
after using the MS renormalization prescription and choosing the renormalization parameter p, so that 
the zero temperature gap equation is satisfied. 


we have the following expression for the pressure (in units of Aq), 
P{T) 


(56) 




= 3 


/(^A'^ r, T') + /(-zA' 2, r, T') - -/(O, r, T') 


3 
2 
9A'4 
327r2 


/(^A'^ 0, T') + 0, T') - -1(0,0, T') 


InA'^ - - - 


647r2 


(57) 


In (57) prime quantities stand for quantities in units of Aq, while A and Aq satisfy their gap equation. 
The last term of (57) accounts for the zero temperature subtraction, so that p(0) = 0, according to the 
definition of I{m‘^,r,T) in (40). Note that the coupling constant does not explicitly appear in (57) and 
that Ao stands for the Gribov parameter at T = 0. 


The interaction measure / is defined as the trace anomaly in units of T^, and I is nothing less than 
the trace of the of the stress-energy tensor, given by 


dfiu = {p + f)u^Ui, - prj^^ , (58) 

with e being the internal energy density, which is defined as e = £y Ts (with s the entropy density), 
u = (1,0,0, 0) and the (Euclidean) metric of the space-time. Given the thermodynamic definitions 
of each quantity (energy, pressure and entropy), we obtain 

I = . (59) 
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Both quantities display a behavior similar to that presented in |39j (but note that in they plot the 
temperature in units of the critical temperature (T^ in their notation), while we use units Ag). Besides 
this, and the fact that we included the effect of Polyakov loop on the Gribov parameter, in [33] a 
lattice-inspired effective coupling was introduced at finite temperature while we used the exact one-loop 
perturbative expression, which is consistent with the order of all the computations made here. 

However, we notice that at temperatures relatively close to our Tc, the pressure becomes negative. 
This is clearly an unphysical feature, possibly related to some missing essential physics. For higher 
temperatures, the situation is fine and the pressure moreover displays a behaviour similar to what is 
seen in lattice simulations for the nonperturbative pressure (see for the SU{3) case). A similar 
problem is present in one of the plots presented in |391 Fig. 4], although no comment is made about 
it. Another strange feature is the oscillating behaviour of both pressure and interaction measure at 
low temperatures. Something similar was already observed in m where a quark model was employed 
with complex conjugate quark mass. It is well-known that the gluon propagator develops two complex 
conjugate masses in Gribov-Zwanziger quantization, see e.g. [HlinilMlITTj for some more details, so 
we confirm the findings of m that, at least at leading order, the thermodynamic quantities develop an 
oscillatory behaviour. We expect this oscillatory behaviour would in principle also be present in |39] if 
the pressure and interaction energy were to be computed at lower temperatures than shown there. In 
any case, the presence of complex masses and their consequences gives us a warning that a certain care 
is needed when using GZ dynamics, also at the level of spectral properties as done in 174||75|, see also 

[ISIIZS]- 


These peculiarities justify giving an outline in the next Section of the behaviours of the pressure and 
interaction measure in an improved formalism, such as in the Refined Gribov-Zwanziger one. 




Figure 4: Left: GZ pressure (relative to the Stefan-Boltzmann limit pressure ~ T^). Right: GZ trace 
anomaly. 


6 Outlook to the Refined Gribov-Zwanziger formalism 


The previous results can be slightly generalized to the case of the Refined Gribov-Zwanziger (RGZ) 
formalism studied in [1111311131731110]. In this refined case, additional nonperturbative vacuum con¬ 
densates such as (A^) and are to be introduced. The corresponding mass dimension two op¬ 

erators get a nonzero vacuum expectation value (thereby further lowering the vacuum energy) and thus 
influence the form of the propagator and effective action computation. The predictions for the RGZ 
propagators, see also miiHiiHi], are in fine agreement with ruling T = 0 lattice data, see e.g. also 
mi m m limij nmol |3TJ m. THs is m contrast with the original GZ predictions, such that it could 
happen that the finite temperature version of RGZ is also better suited to describe the phase transition 
and/or thermodynamical properties of the pure gauge theory. 

Due to the more complex nature of the RGZ effective action (more vacuum condensates), we will 
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relegate a detailed (variational) analysis of their finite temperature counterpart^ to future work, as this 
will require new tools. Here, we only wish to present a first estimate of the deconfinement critical tem¬ 
perature Tc using as input the T = 0 RGZ gluon propagator where the nonperturbative mass parameters 
are fitted to lattice data for the same propagator. More precisely, we use m 

Aab( . cab p^+M'^+pi / 

(M2 + + pY) + m2 (#2 -tpi )■+ A4 \ p2 

where we omitted the global normalization factor Z which drops out from our leading order computatior0 
In this expression, we have that 









b ,„ab 
p, Vp ■ 


-ab-ab 

■ ^p ^p 


) Pi 


(61) 


The free energy associated to the RGZ framework can be obtained by following the same steps as in 
section leading to 


£viT) 


(d-1) 




(d-1) 


-f 


d'^p 




In 


T) + I{rl,r, T) - I{N\ r, T) - r, T) 

a — 1 

0, T) + /(r^ , 0, T) - /(TV", 0, T) - 0, T) 

p* + (m" -I- N'^)p^ + + A^) A SA^d 

p2 -(- ]V2 J 45'2 ’ 


(62) 


with standing for minus the roots of the denominator of the gluon propagator (601, iV" = M" + pi , 
and I{ni^,r,T) given by (40). Explicitly, the roots are 


2 (m" -I- N"^) ± a/( m2 -I- A^ 2)2 _ 4 (j^ 2 jY 2 _|_ 


The (central) condensate values were extracted from [77] : 


TV" = M"-fpi =2.51GeV" , 
m" = -1.92 GeV" , 
A'‘ = 5.3GeV^ . 


(63) 


(64a) 

(64b) 

(64c) 


Once again the vacuum energy will be minimized with respect to the Polyakov loop expectation value 
r. For the analysis of thermodynamic quantities, only contributions coming from terms proportional to 
I{jm?,r,T) will be needed. Therefore, we will always consider the difference £v{T) — £y{T = 0). Since 
in the present (RGZ) prescription the condensates are given by the zero temperature lattice results (64) 
instead of satisfying gap equations, the divergent contributions to the free energy are subtracted, and no 
specific choice of renormalization scheme is needed. Furthermore, explicit dependence on the coupling 
constant seems to drop out of the one-loop expression, such that no renormalization scale has to be chosen. 
Following the steps taken in Section [5.1| we find a second order phase transition at the temperature: 


T,,Ht = 0.25GeV 


(65) 


which is not that far from the value of the SU(2) deconfinement temperature found on the lattice: 
Tc « 0.295GeV, as quoted in [Ml 157] . 

In future work, it would in particular be interesting to find out whether —upon using the RGZ 
formalism— the Gribov mass and/or RGZ condensates directly feel the deconfinement transition, similar 
to the cusp we discovered in the Gribov parameter following the exploratory restricted analysis of this 
paper. This might also allow to shed further light on the ongoing discussion of whether the deconfinement 
transition should be felt at the level of the correlation functions, in particular the electric screening mass 
associated with the longitudinal gluon propagator [M] IMJ 110011101] . 

^From | 93l 194119^ . the nontrivial response of the d = 2 condensate {A^) to temperature already became clear. 

®This Z is related to the choice of a MOM renormalization scheme, the kind of scheme that can also be applied to lattice 
Green functions, in contrast with the MS scheme. 
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Let us also consider the pressure and interaction measure once more. The results are shown in Figure 
and Figure respectively. The oscillating behaviour at low temperature persists at leading order, 
while a small region of negative pressure is still present — see the right plot of Figure These findings 
are similar to |102] (low temperature results are not shown there), where two sets of finite temperature 
RGZ fits to the SU{3) lattice data were used |10dl ll()4| . in contrast with our usage of zero temperature 
SU{2) data. In any case, a more involved analysis of the RGZ finite temperature dynamics is needed to 
make firmer statements. As already mentioned before, there is also the possibility that important low 
temperature physics is missing, as for instance the proposal of |102j related to the possible effect of light 
electric glueballs near the deconfinement phase transition |105L I106j . Obviously, these effects are absent 
in the current treatment (or in most other treatments in fact). 



Figure 5: Right and left plots refer to the RGZ pressure in terms of T/Tc and in units of T^. In the left 
plot, a wide temperature range of is shown. In the right plot, a zoom is made for temperatures around 
1.10 Tc to show the existence of negative pressure. 


l/T‘* 



Figure 6: The RGZ interaction measure I/T^ in units T/Tc- 


7 Summary 

In this paper we studied the Gribov-Zwanziger (GZ) action for SU{2) gauge theories with the Polyakov 
loop coupled to it via the background field formalism. Doing so, we were able to compute in a simultaneous 
fashion the finite temperature value of the Polyakov loop and Gribov mass to the leading one-loop 
approximation. The latter dynamical scale enters the theory as a result of the restriction of the domain 
of gauge field integration to avoid (infinitesimally connected) Gribov copies. Our main result is that 
we found clear evidence of a second order deconfinement phase transition at finite temperature, an 
occurrence accompanied by a cusp in the Gribov mass, which thus directly feels the transition. It is 
perhaps worthwhile to stress here that at temperatures above Tc, the Gribov mass is nonzero, indicating 
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that the gluon propagator still violates positivity and as such it rather describes a quasi- than a “free” 
observable particle, see also pSl 1107] for more on this. 

We also presented the pressure and trace anomaly, indicating there is a problem at temperatures 
around the critical value when using the original GZ formulation. We ended with a first look at the 
changes a full-fledged analysis with the Refined Gribov-Zwanziger (RGZ) formalism might afflict, given 
that the latter provides an adequate description of zero temperature gauge dynamics, in contrast to the 
GZ predictions. This will be studied further in upcoming work. Note that, even not considering finite 
temperature corrections to the condensates in the RGZ formalism, the region of negative pressure is 
considerably smaller than the region found with the GZ formalism. 

A further result of the present paper, which is interesting from the methodological point of view, is 
that it shows explicitly that finite-temperature computations (such as the computation of the vacuum 
expectation value of the Polyakov loop) are very suitable to be analyzed using analytical Gasimir-like 
techniques. The interesting issue of Gasimir-style computations at finite temperatures is that, although 
they can be more involved, they provide one with easy tools to analyze the high and low temperature limits 
as well as the small mass limit. Moreover, within the Gasimir framework, the regularization procedures 
are often quite transparent. Indeed, in the present paper, we have shown that the computation of the 
vacuum expectation value of the Polyakov loop is very similar to the computation of the Gasimir energy 
between two plates. We believe that this point of view can be useful in different contexts as well. 
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